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Comment on "Density and Spin response of a strongly- 
interacting Fermi gas in the attractive and quasi-repulsive 
regime" 

In Ref.jjTl] the authors summarize a linear response the- 
ory for the superfluid Fermi gases undergoing BCS-BEC 
crossover. They include a rather complex set of diagrams on 
the basis of avoiding a divergence in the density response. We 
wish to point out that (1) this divergence is associated with an 
unphysical first order transition in their theory (2) repairing 
this divergence, as they do, leads to an unphysical Meissner 
effect in the normal state and therefore a violation of the well 
known transverse f-sum rule. (3) In earlier work we have pre- 
sented in a series of papers yUg] (not cited by the authors), 
a systematic theory of spin and charge linear response above 
and below T c in a strongly-interacting (attractive) Fermi gas 
which does not suffer from concerns (1) and (2). In addition 
(4) the material presented in Reference [ 1 ] and the Supple- 
ment suggests that it is relatively straightforward to extend the 
authors' calculations below T c . We assert that this statement 
is misleading in the density channel, because of the necessity 
of including collective mode effects, not discussed by the au- 
thors. 

We confine our attention in this Comment to the attractive 
interaction regime which was discussed by us in earlier work 
We first show that there will be a Meissner effect in 
the normal state of Ref. A consistent set of diagrams for 
linear response must necessarily be related to the self energy 
£ through the Ward Identities. In Appendix D of Ref. [2] 
we have discussed the appropriate set of transport diagrams 
associated with the self energy used in the Nozieres Schmitt- 
Rink (NSR) approach and in Ref. lH. 

The superfluid density is given by — = ( — ) — 

r J D J m Km J xx 

Pxx{0) where P X x(0) denotes the transverse component of 
the current-current response function in the limit of zero fre- 
quency and zero momentum. A properly consistent theory 1 2] 
for this correlation function involves the fully dressed Green's 
function G{K), where if is a 4-vector, the free fermion dis- 
persion £fc and diagrams of the Maki-Thompson (MT) and 
Aslamazov Larkin (AL) type. We refer to the bare vertex as A 
associated with the density of states (DOS) diagram. In order 
for there to be no unphysical Meissner effect above T c , one 
must take 


This additional set of diagrams, argued to avoid 
an unphysical divergence in compressibility, leads 
in the normal phase to a finite superfluid density 
n s /m = -2J2 K G 2 {K)AL new (K,K) x \{K,K) x , which 
results from the inconsistent approximations in the self 
energy and the response functions. 

While the authors of Ref. J3] suggest that it is straight- 
forward to extend their compressibility calculations below 
the transition temperature, this statement is incorrect. In the 
broken-symmetry phase it has long been known that a consis- 
tent calculation of the compressibility ^ based on the lon- 
gitudinal density-density correlation functions will require, 
through gauge invariance, that the collective modes of the or- 
der parameter are also included [3]. Importantly, they are very 
complicated to address away from the BCS regime. 

In summary, a failure to maintain gauge invariance leads to 
the presence of a Meissner effect in the normal state. Since we 
have just noted that the compressibility obtained via linear re- 
sponse theory is extremely sensitive to gauge invariance, one 
should be concerned about the validity of any theory of the 
compressibility which is associated with a normal state Meiss- 
ner effect. Interestingly enough, the normal state Meissner 
effect in Ref. yj] is a consequence of the authors' attempt to 
regularize the compressibility. 

This Comment is intended to emphasize the key role played 
by sum rules and other conservation constraints in transport 
and scattering theories of the Fermi gas superfluids, as we 
have repeatedly stressed in our (uncited) papers [2-6;]. The 

ood quantitative agreement with experiment reported in Ref. 

l|] is not a valid reason for ignoring theoretical inconsisten- 
cies. 
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P£f H (0) = —2 G (K)[X(K, K) + MT(K, K) 
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While these were all included in Ref. [1], additional diagrams 
of the AL type also appeared (Fig. 1(e) of Ref. jlJ])- Let 
AL new denote this set of diagrams. Then the current-current 
response function from the diagrams selected in Ref. |l|] is 

P a p(0) 

= I% SR {0) - 2 ° 2 ( R ) AL ne W (K, K) a X(K, K) p . (2) 
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